Nonsingular vortices in (s+d)-wave superconductors 
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The structure of a single flux line in (s+d)-wave superconductors has been analyzed within the Ginzburg- 
Landau (GL) model generalized for two order parameter components. The fourfold symmetric singular vortex 
solution is shown to be unstable in a certain range of the GL parameters with respect to the mutual shift of s- and 
d- wave unit vortices. The resulting nonsingular vortex structure is studied both analytically and numerically. 



Recently the distinctive characteristics of vor- 
tices in unconventional superconductors which 
can be described by the phenomenological 
Ginzburg-Landau (GL) theory with a multicom- 
ponent order parameter (OP) are of great in- 
terest in connection with the investigations of 
the mixed state structure in high-T c and heavy 
fermion compounds which are strong candidates 
for unconventional superconductivity. A single 
flux line in such systems is known to contain a 
set of unit vortices of different OP components 

Q . One can indentify two possible types of 
flux lines: (i) singular vortices (which have at 
least one point where the superconducting gap 
is zero) and (ii) nonsingular vortices (where the 
gap is nonzero everywhere in the vortex core). 
In this paper we focus on the case of s + d x 2_ y 2- 
wave pairing (which can be relevant to the case of 
high-T c superconductors) and consider the range 
of GL functional parameters where the singular 
vortices (studied in j^-Q]) become unstable ac- 
cording to the scenario analogous to the one pro- 
posed in 0J^] for heavy fermion compounds. The 
goal of this paper is to analyse the detailed struc- 
ture of resulting nonsingular vortices using both 
numerical and analytical methods. 

We start with the GL free energy functional 
generalized for two components of the OP ^> d and 
ty s corresponding to the dj^^-wave and s-wave 
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pairing, respectively 
F : / {a, >!'/-' ■ a s >l'. 2 

(3, 



b 4\*,U b ^ 



-a\v d \ 2 \* s \ 2 +^ 2 d *f+*f* 2 s ) 
-K s \n^ s \ 2 +K d \n^ d \ 2 + 7 [(u* x ^*u x ^ d 



U* y ^* s Uy tfrf) + C.C.] 



— }dr, 



(1) 



where II = V - i^A, H = curl A, r = (x,y), 
a s =a s (T - T cs ), a d = a d (T - T cd ) (we assume 
T cs < T cd and the magnetic field is applied paral- 
lel to the c-axis). If the GL parameters are chosen 
so that there is an one-component homogeneous 
state (H = 0), then for H ^ the subdominant 
OP can appear in the vortex core regions with the 
inhomogeneous dominant OP component either 
due to the mixing gradient terms (7 =/= 0) [^] or 
due to the instability of the inhomogeneous state 
against the formation of the subdominant OP nu- 
cleus Jl],|],||. Let us consider the latter mecha- 
nism (which is obviously responsible for the for- 
mation of nonsingular vortices) for the simplest 
case 7 = 0, (3 = and continue with the analysis 
of the linearised GL equation of \I/ S : 



K s V 2 * s +a\* d \ 2 * s 



(2) 



where ty d describes the vortex solution within 
the conventional single-component GL theory, 
also we neglect A. The "lowest energy eigen- 
state" of the Schrodinger equation (0) defines 



the temperature T* of the phase transition into 
nonsingular vortex state. Using the following 
approxima tion = \$> d\ cor / \/ r 2 + 2^ d 2 (here 
|*dU = \f\aJJb~d, £,d = y/K d /\a d \) for the cases 
K s <C K d a s /a d and a <C bdK s /Kd, we obtain 
-a(T*) = y/2K s a/(K d b d ) and -a(T*) ~ a/6 d) 
respectively (here a(T) = a s (T)/|a d (T)|). We 
propose the following extrapolation for the phase 
transition curve on a — T plane (Fig.|l|), which 
separates the regions with singular and nonsingu- 
lar vortices: a = -b d a{T*){l - a{T*)K d /{2K s )). 
This expression is valid for the above two limiting 
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Figure 1. The phase diagram on a — T-plane 
for a s = 3.2a d , T cs = 0.9T crf , b s = % d , = 0, 
K s — Q.SKd- The homogeneous states with 
* s = 0, *d ^ and with ^> d ^ exist above 
and below the dashed line, respectively. The non- 
singular vortices exist between the dashed line 
and the solid line. 



cases and for the particular case (a s = 1.25ad, 
T = 0, T cs = 0.8T cd , b s = b d , 13 = 0, K s = K d , 
7 = 0) studied in (l). The analysis of the struc- 
ture of nonsingular vortices has been carried out 
using numerical calculations based on the time- 
dependent GL theory. It was found out that 
for 7^0 the nonsingular vortices with broken 
fourfold symmetry (see Fig.||) are stable below 
the phase transition curve (solid line in FigJ^), 
whereas above the curve it is the singular four- 
fold symmetric vortices that are energetically 
favourable. 



Figure 2. The typical structure of s-wave OP for 
nonsingular vortex for a s = 3.2a d , T cs = 0.9bT cdl 
T = 0.68T cd , b s = 9b d , a = 3b d , (3 = 0, 
K s = 0.5K d , 7 = 0.65K d . The inset shows the 
positions of d-wave vortex whose winding num- 
ber N = +1 (black diamond), the s-wave vortices 
with N = +1 (circles) and with N = — 1 (cross). 



In conclusion, we have considered the structure 
of nonsingular vortices in s + d x 2_ y 2 supercon- 
ductor and found the range of parameters where 
these vortices exist. The nonsingular vortices in 
such materials have neither the fourfold symme- 
try, nor the normal core region. This can lead to 
some interesting effects (e.g. nontrivial dynamics, 
moderate pinning effects etc.). 
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